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We describe the effect of power-law initial distributions of clusters on ordinary percolation and its 
generalizations, specifically, models of explosive percolation processes based on local optimization. 

These aggregation processes were shown to exhibit continuous phase transitions if the evolution 
starts from a set of disconnected nodes. Since the critical exponents of the order parameter in 
explosive percolation transitions turned out to be very small, these transitions were first believed 
to be discontinuous. In this article we analyze the evolution starting from clusters of nodes whose 
sizes are distributed according to a power law. We show that these initial distributions change 
dramatically the position and order of the phase transitions in these problems. We find a particular 
initial power-law distribution producing a peculiar effect on explosive percolation, namely before 
the emergence of the percolation cluster, the system is in a “critical phase” with an infinite general¬ 
ized susceptibility. This critical phase is absent in ordinary percolation models with any power-law 
initial conditions. The transition from the critical phase is an infinite order phase transition, which 
resembles the scenario of the Berezinskii-Kosterlitz-Thouless phase transition. We obtain the crit¬ 
ical singularity of susceptibility at this peculiar infinite-order transition in explosive percolation. It 
turns out that susceptibility in this situation does not obey the Curie-Weiss law. 

PACS numbers: 64.60.ah, 05.40.-a, 64.60.F- 


I. INTRODUCTION 

The percolation transition is the key phase transition 
occurring in disordered systems including disordered lat¬ 
tices and random networks EHi. The gradual increase 
of the number of links in a network or lattice leads to 
the growth of clusters of connected nodes and eventu¬ 
ally to the formation of a percolation cluster (giant con¬ 
nected component) at the percolation threshold. This 
phase transition was studied in detail, and understood 
to be continuous in all disordered systems which were 
explored. Recently a new class of irreversible percolation 
processes, so-called “explosive percolation”, was intro¬ 
duced 0, where the new links are added to the sys¬ 
tem using metropolis-like algorithms. Although these 
processes directly generalize ordinary percolation, they 
demonstrate a set of features remarkably distinct from or¬ 
dinary percolation. The unusual properties of this kind 
of percolation led to the initial reports based on sim¬ 
ulations [Ml that these processes show discontinuous 
transitions. Solving the problem analytically, we have 
shown that the explosive percolation transitions are ac¬ 
tually continuous [Q . These transitions have a so small 
critical exponent of the percolation cluster size, that in 
simulations of finite systems they can be easily perceived 
as discontinuous [IJ, [l^ . This conclusion was supported 
by subsequent works of physicists jl^-[l^ and mathe¬ 
maticians [^ . 

In our previous work [2lj we developed the scaling the¬ 
ory of explosive percolation phase transitions for a wide 
range of models, explaining the continuous nature of the 
transitions and their unusual features. We obtained the 
full set of relevant critical exponents and scaling functions 


in the typical situation, in which the evolution starts from 
isolated nodes or clusters with sufficiently rapidly decay¬ 
ing size distribution. In our papers [14l.ll5l.l2ll| we employ 
the following model of explosive percolation, at which the 
number of nodes N is fixed and links are added one by 
one. At each step we choose two sets of m random nodes, 
from each set we select the node that is in the smallest 
of m clusters, and then we add a new link between these 
two nodes. 

In the present article we show that the explosive perco¬ 
lation transition, as well as ordinary percolation, strongly 
depends on the initial conditions of the process. In par¬ 
ticular, slowly decaying initial cluster size distributions 
can change crucially the nature of these transitions. The 
effects are interesting and add much to understanding of 
explosive percolation and other generalizations of ordi¬ 
nary percolation. So in the present article we explore in 
detail the effect of initial conditions on the percolation 
transitions in systems including ordinary and explosive 
percolation models. We consider power-law initial clus¬ 
ter size distributions with exponent f, and for different 
values of the exponent find a spectrum of distinct criti¬ 
cal behaviors. Here we introduce the initial cluster size 
distribution exponent f in contrast to the critical clus¬ 
ter size distribution exponent traditionally denoted by r. 
Because of the power-law critical distribution we expect 
that power-law initial conditions produce interesting ef¬ 
fects. We will indicate the range of f where the transition 
point coincides with the initial moment of the process, 
tc = 0. In particular, for ordinary percolation, if f = 3, 
the percolation cluster emerges as S' ^ exp(—const/t), 
where S is the relative size of this cluster. In contrast, 
for explosive percolation, we find that there exists a value 
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f < 2-t 1/(2 to- 1) 

f = 2-tl/(2TO-l) 

T > 2-t 1/(2 to - 1) 


tc = 0 

tc = 0 

tc > 0 


P(s,0) 

P(s,0) ~ 

P(s,tc) ~ 

to = 1 


S ~ exp (-const/t) 

S^(t-tcf 


P = {t- 2)/(3 - t) 

— 

/3 = 1 



x~t-" 

X ~ |t - tc|“^ 


tc = 0 

tc > 0 

tc > 0 


P(s,0) 

P(s, tc) ~ 3^“” In^ s 

P{s,tc) ~ 



S ~ exp[-const/(t — tc)'^] 

S ~ (t - tc)'’ 

TO > 1 

/3= 

^ 1-(2to- !)(-?- 2) 

1 

^ ^ A(2to- 1) - 1 

r - 2 ~ ^ ~ g-1.43m ^ 


x~t-i 

X ~ (t — tc)“^“'^ if t > tc 

X ~ |t - tcT’ 



X = oo if t < tc 



TABLE I. Summary of results. The initial distribution P{s,t = 0) ~ f > 2, 5 is the relative size of the percolation 

cluster, and x is the susceptibility. At r = 2 + l/(2m — 1), we show only the most singular factor of S. 


of f at which the phase transition turns out to be infinite- 
order and occurs at tc > 0. In this situation, the system 
a,t t < tc is in the “critical phase” with divergent sus¬ 
ceptibility. We also find susceptibility at t > tc for any 
m and show that its critical exponent is nontraditional, 
differing from the Curie-Weiss law. The main results of 
the paper are presented in Table HI 

The paper is organized in the following way. Section ITTl 
outlines our results and methods. In Sec. uni we consider 
effect of power-law initial conditions on ordinary perco¬ 
lation (to = 1), which is the simplest particular case of 
the more general model analyzed in this work. In Sec. lIVI 
we study the effect of initial conditions on the explosive 
percolation model (to > 2), which turns out to be prin¬ 
cipally different from the case of to = 1. 


II. RESULTS 

To help the reader let us outline the main results of 
this article. We use the following set of models. At each 
time step a new link connecting two nodes is added to the 
network of N nodes. At each step sample two times: (i) 
choose TO > 1 nodes uniformly at random and compare 
the clusters to which these nodes belong; select the node 
within the smallest of these clusters; (ii) similarly choose 
the second sampling of to nodes and, again, as in (i), 
select the node belonging to the smallest of the to clus¬ 
ters; (iii) add a link between the two selected nodes thus 
merging the two smallest clusters. The resulting process 
is described by the time dependent probability P{s,t) 
that a randomly chosen node belongs to a finite cluster 
of size s, where the time t = L/N, where L is the num¬ 
ber of added links (number of steps of the process). We 
assume that N is infinite. Then this aggregation process 


is described by the evolution equation 

—^^=s^Q{u,t)Q{s-u,t)-2sQ{s,t), (1) 

U — 1 

where Q{s,t) is the probability that a cluster chosen to 
merge is of size s. This probability is expressed in terms 
of P{s,t). In particular, when to = 1, the distribution 
Q{s,t) coincides with P{s,t), and the model is reduced 
to ordinary percolation. 

We focus on the effect of slowly decaying initial distri¬ 
butions P(s, t = 0), namely power laws P{s, 0) = aos^“”, 
which can produce tc = 0 or, in the case of explo¬ 
sive percolation, a critical phase. Using the evolution 
equation © we analyze the Taylor expansion P(s,t) = 
P{s, 0) -I- Ai(s)t -I- A 2 t^ -I- ... and obtain the scaling form 
of the distribution P{s, t) = s^“'”/(st^/‘^). We obtain the 
relation between the scaling function f{x) and the size 
of the percolation cluster S ^ . We demonstrate that 

when the exponent /3 is integer, there are two contribu¬ 
tions to S. The first is the well known contribution deter¬ 
mined by the scaling part of P(s, t) (see, for example, the 
book i)- However, there is a second, analytic, contri¬ 
bution that was not considered in Ref. 0- In particular, 
the exponent /3 is 1 for ordinary percolation with rapidly 
decaying P(s, 0). In this standard case, it is the combina¬ 
tion of these two contributions S{t) = const (|t—tc|+t—tc) 
that produces S{t < tc) = 0 and the proper dependence 
S{t > tc). 

In Table U we present a summary of the main results 
of this article. The first row of the table shows the criti¬ 
cal behaviors of the ordinary percolation model, to = 1, 
for different f > 2. The values f > 3 result in a perco¬ 
lation phase transition at tc > 0 with standard critical 
exponents. On the other hand, at 2 < f < 3 the per¬ 
colation cluster emerges at t = 0 with critical exponents 
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given in terms of f. In the marginal case of f = 3 we ob¬ 
serve an infinite-order phase transition with a singularity 
S ^ exp(—const/t). This set of the critical singularities 
of S' for m = 1 agree with those obtained in Ref. [2^ . 
Interestingly, the susceptibility y = (s)p = 
obeys the Curie-Weiss law y ^ {t — tc)~^ in all the con¬ 
sidered situations except at f = 3, when y ~ (note 
that f < 3 gives tc = 0). 

The second row of the table presents our results for ex¬ 
plosive percolation m > 1. For these m, an infinite order 
phase transition occurs when r = 2-|-l/(2m—1) at tc > 0. 
Before tc the system is in the critical phase, in which 
the susceptibility [2l| diverges, while in the percolation 
phase the critical behavior of susceptibility y(< > tc) dif¬ 
fers from the Curie-Weiss law. We find that the size 
distribution of clusters at the point of the infinite-order 
phase transition is P{s,tc) s^“'^ln^s. We obtain the 
exponent A close to 1 solving the evolution equation nu¬ 
merically. 

Finally, we obtain a general relation between the sus¬ 
ceptibility and the size of the percolation cluster close to 
the critical point. 


mdlnS 
T dt 


( 2 ) 


This simple relation is valid for all models and initial 
conditions considered in this work. The detailed deriva¬ 
tions of these analytical results are given in subsequent 
sections. 


III. EFFECT OF INITIAL CONDITIONS IN 
ORDINARY PERCOLATION (m = 1) 


whose solution p{z,t) can be obtained from 

lnz = 2t{l- p)+g{p), (6) 

where the function g{p) is determined by initial condi¬ 
tions. 

Let us find this function. Our results will be completely 
determined by the asymptotic of the initial distribution. 
An initial cluster size distribution with a power-law tail, 
i.e., P{s) = aos^~'^ for large s, leads to the following 
singularity of the generating function 

OO 

p{l) - p{z)= ^ P(s)(l - = ao E 

S—1 S 

nCO 

= ao(—lnz)'^“^ / dx x^~'^{1 — e~^) 

Jo 

=-aor(2 - r) (-lnz)^“^ , (7) 


at z = 1. 

We find the function g{p), replacing the left-hand side 
of Eq. ([6]) by Eq. ([T]) and putting p{l) = 1, since S = 
1 — p(l) = 0 at t = 0. Then Eq. ([6]) becomes 


which is valid when z —>■ 1. Putting z = 1 and 1 — p = S 
we get 


2tS = 


^ ) 

aor(2-T)J 


l/(f-2) 


(9) 


In the case of m = 1, our process is actually ordi¬ 
nary percolation, that is at each step we chose uniformly 
at random two nodes and interconnect them. This can 
be treated as an aggregation process in which at each 
step two clusters, chosen with probability proportional 
to their sizes, merge together. This process is described 
by the probability distribution P(s,t) that a randomly 
selected node belongs to a cluster of size s at moment 
t (each step increases time t by 1/iV) . In the infinite 
system, the evolution of this distribution is described by 
the following Smoluchowski equation [2^ : 

= s^P{u,t)P{s-u,t)-2sP{s,t) (3) 

U—1 

for a given initial distribution P{s,t = 0). Defining the 
generating function p{z,t) as 

OO 

P{z,t) = '^P{s,t)z\ (4) 

S = 1 


we can rewrite Eq. ([3]) in terms of p, and analyze the 
resulting partial differential equation: 


dp{z,t) 


2[l-p(z,t)] 


dp{z,t) 


(5) 


The last equation has two solutions, the trivial one S = 0 
and a nontrivial solution 

S = (2t)(^-2)/(3-D [-aor(2 - ^ (10) 

If f > 3 then r(2 — f) > 0, and Eq. (jH]) has only one 
real solution S = 0, showing that the transition does not 
occur at tc = 0 for this range of f. In this case, to find a 
real solution S' > 0 for t > tc > 0 we must also consider 
the analytic terms omitted in Eq. ©. Eor the range 
2 < r < 3, we have r(2 — t) < 0 and the solution (flOl) is 
real and positive, that is 

^^(f-2)/(3-f)^ (II) 

for small t, and B — [— aor(2 — 

Therefore, if r < 3, the transition occurs at the initial 
moment. Note that for r < 3, the first moment of the ini¬ 
tial distribution diverges. For ordinary percolation this 
moment has the meaning of susceptibility Q, and its di¬ 
vergence indicated that the poit t = 0 is indeed the crit¬ 
ical point. We will consider the case of f = 3 separately 
in Sec. IIII Bl 

It is easy to see that the distribution P{s, t) is an ana¬ 
lytic function of t. (i) The initial distribution P{s, 0) has 


dt 


d\nz 
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no divergencies at any s. (ii) Let us, for a moment as¬ 
sume that P{s,t) at has singularity with non-integer 
(j) at t = 0. Then the lowest non-integer power on the 
left-hand side of Eq. m is while on the right-hand 
side the lowest non-integer power is which shows that 
the assumption was not correct. Then we can write the 
Taylor expansion of the function P{s,t) around t = 0: 

P{s,t) = Aq (s) -|- Ai (s)t -|- A 2 A 3 . (12) 


The first term in expansion m is the initial distribu¬ 
tion, Ao(s) = P(s, 0). The coefficient of the second term, 
Ai(s), is the first derivative dtP{s,t = 0 ), ^ 2 ( 5 ) is the 
second derivative {l/2)d^P{s,t = 0), and so on. Then, 
given an initial distribution P(s, 0), we can find the co¬ 
efficients Ai(s) sequentially differentiating both sides of 
Eq. @. We analyze these coefficients in different ranges 
of f. 

The remainder of this section is organized in the fol¬ 
lowing way. In Sec. IIII AI we consider the case of f < 3. 
In this region, we derive the scaling of the distribution 
P(s, t) containing a scaling function. We obtain the rela¬ 
tion between the analytical features of this scaling func¬ 
tion and the singularity of the relative size S of the perco¬ 
lation cluster at t = 0. In Sec. IIII Bl we consider the case 
of f = 3. We show that in this situation all the deriva¬ 
tives of S{t) with respect to t are zero at t = 0 and derive 
the respective scaling function. In Sec. IIII Cl we analyze 
the singularity of the susceptibility of this problem, and 
its relation with the percolation cluster size S. 


A. The case of f < 3 

Let us first consider the case off < 3, derive the scaling 
of the distribution P{s,t), and the critical singularity of 
the percolation cluster size. 


1. Scaling of P{s,t) 


We consider initial configurations without a percola¬ 
tion cluster, so X]s^o(s) = I- We find the coefficients 
A„(s), n> 0 , using the identity dl'^^P{s, t = 0) = A„(s) n! 
and Eq. ©. For the second coefficient we have 

s-l 

Ai(s) = s ^ Ao('u)Ao(s - u) - 2sAo(s). (13) 

U—1 


The sum X]jj<s ^o('a)Ao(s — u) cannot be directly re¬ 
duced to the integral dx [a:(I —for large 

s, because it diverges at both limits 0 and I. So, following 


our work we rewrite Eq. m as 

s-l 

^i(s) [Ao(u)-Ao(s)] [Ao(s-u)-Ao(s)] 

U—1 

S-l 

- s(s-I)[Ao(s)]^ -I- 2Ao(s) ^ Ao{u) - 2sAo(s) 


S-l 

["4o(m)-Ao(s)] [Ao(s-u)-Ao(s)] 

U—1 


00 

- s(s-I)[Ao(s)]^ - 2sAo(s) ^ Ao(u), (14) 

U—S 


where we took into account the normalization condition 
Ao(s) = I. For large s, the sums in last equation can 
be already reduced to integrals, which converge at both 
limits, and Ao(s) can be replaced by Then we 

get 

Ai(s) = f dx \x^~'^ — l] [(I — xY~'^ — l] 

Jo 

/ OO 

dxx^-^ 


= a 


2 

0 


r(4 - 2f) 


= ais 


4-2f 


(15) 


where we have introduced the coefficient ai. 

The case off = 5/2 is special. For this f the coefficient 
oi in Eq. (flSj) is zero, and Ai decays faster than To 

find the asymptotics of Ai in this situation, we must take 
into account the higher-order terms that were neglected 
when passing from the sums in Eq. da to the integrals 
of Eq. (fT31) . In general, we estimate the difference of the 
respective integral and sum with arbitrary exponents tp 
and (p 


s du[u ^—s ’^][(s—u) s 
Jo 


s-l 

-s^[u-’^-s-’^][(s-u) 

u—1 


-s-'f’] = 


(16) 


where we assume that Q < ip^tp < 2. In particular, for 
Ai, we have ip = (p = f — 1. Then, when f = 5/2, the 
large s asymptotics of Ai is 

Ai(s) oc —oc —if f = 5/2. (17) 

In Appendix 1X1 we calculate the asymptotics of A 2 (s): 

For the particular values f = 7/3 and 8/3 the coefficient 
02 = 0. To find the asymptotic behavior of A 2 in these 
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cases it is necessary to consider higher-order terms that 
were neglected when passing from Eq. (IA2I) to Eq. (IMI) . 
similarly to how we treated Ai. This analysis gives 


^ 2 ( 5 ) oc 


— ^ (X —s 

— OC — 


if f = 7/3, 
if T = 8/3. 


(19) 


These calculations can be repeated for the next terms 
of expansion and in general we find that the n-th 
term 


^n('5) 


= a„s 


1 —f+n(3—f) 


( 20 ) 


for large 5 . In Appendix [B1 we obtain the general expres¬ 
sion 


[2aor(2-f)]"+^ 

2r[(n-fl)(2-f)](n-fl)!' 


( 21 ) 


Notice that this expression generalizes the results for oi 
and 02 , Eqs. m and (I18|) . respectively. The coefficients 
a„ are expressed in terms of oq and f and become zero 
for f = 3 — l/(n-|-l). For this f, the coefficient An,(s) 
decays as 

Equation (|20|) enables us to write the Taylor expansion 
of P{s,t) in the form; 


P{sA) = s^ = A), (22) 

n=0 


The function f(x) is the scaling function of the prob¬ 
lem. This function is only analytic at zero for t = 'bj^ 
when coefficients a„ are zero for odd n. In general, for 
2 < r < 3, f{x) is represented as the series 

/(a;) =00-1- oia;'^ -I- -I- oax^'^..., (23) 

where 



FIG. 1. Scaling function /(x) calculated using the first 1000 
terms of the series ([231) for oo = 1/and different values of 
f. In inset, /(x) is shown as a function of where 

(T = 3 — f. In this representation the function / is analytic 
at the origin, where its derivative at the origin changes signs 
when r crosses 5/2. 


2. Relation between the singularity of S and the scaling 
function f{x) 

The relative size of the percolation cluster S near tc is 
independent of the details of distribution P{s,tc) in the 
region of small s (non-scaling region). This enables us to 
use the scaling form of distribution P{s,t), Eq. (1^ . to 
recover expression (ED for S{t) near tc- To this end we 
start from the definition of S 

5 = l-^P(s,t), 


cr = 3 — r. 


(24) 


Let us estimate the radius of convergence Xrc of this se¬ 
ries. Expansion (ED is a convergent series for t < r. 


r = lim 

n—¥co 


Anis) 


gl-f-|-ri(3—f) 

^n-l-l(s) 


gl —f-|-(n-|-l)(3—f) 


s^-^. (25) 


Then series (i23l) is convergent for = x < Xrc = 

Let us consider, for example, the particular case r = 
5/2. In this case the scaling function, after substituting 
Eq. (1^ into Eq. (ED: reproduces the known result for 
the mean-field percolation transition at tc > 0 [l|: 




f{x) = gp / ^ 


(—47raQxy 


—ATTaix 

= ape 0 . 


n—0 


(26) 


Note that when f = 5/2 the coefficients a„ in se¬ 
ries of Eq. (05)) with odd n are zero. If we set op 
to in Eq. (12611 . we arrive at the known form 

/(x) = exp(—2x)/-\/^ for the percolation process start¬ 
ing from isolated nodes [l|. Figure [1] shows scaling func¬ 
tions for several values of r. 


— 1 — 'y ( (Ap(s) -|- yli(s)t -|- A2{s)t^ + ...) . (27) 

S 

At t = 0 the relative size of the percolation cluster S' = 0, 
which implies X)s ^o(s) = 1: and so we can write 

s = -^rp„, (28) 

7T,>1 

which is the Taylor expansion of S at zero with the coef¬ 
ficients 


P„ = ^A„(s). (29) 

S 

In Eq. (051) . the coefficients oc They 

diverge ifr < 3— I/(n-|-I) and converge if r > 
3 — I/(u -I- 1). (Recall that at f = 3—l/(n-|-l) the coef¬ 
ficient a„ becomes zero and asymptotics of A„(s) decays 
as For any r < 3 the coefficients Bn are 

infinite for n > ff, where 

/3 = ^- (30) 

3 — r 
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The divergence of these coefficients in the Taylor se¬ 
ries (1^ indicates the singularity of 5 at t = 0. We 
extract this singularity from series (l28)) in the follow¬ 
ing way. We divide the sum X]Tt>i ™to two parts, 
+ Here 

n* = L/3J (31) 

denotes the largest integer smaller or equal to /3. In 
the second part we replace An by their asymptotic form, 
namely 

5 = - ^ Bn^ - E 

l<n<n* n>n* s 

n>n* s 

= - E E 

l<n<n* s n>n* 

-E®”' E o[{st^/nn 

s n>n* 

-- E Bnt--Y.s^-^risf^n- (32) 

l<n<n* s 


all terms leading to divergence, see Eq. (l33l) . Recall¬ 
ing that a = 3 — f, we arrive at the same singularity 
S = 5i(?-2)/(3-f) as in Eqs. (HH), 

B = 2(^-2)/(3-?)[_agr(2 - f)]i/(3-f) 

dxx^~'^ f*{x) — Bn*- (36) 

The first term on the right-hand side of Eq. (l35l) is the 
singular contribution from the scaling behavior P(s, t) = 
s ^~'^The second term on the right-hand side 
of Eq. (1551) is nonzero only when /3 = (f — 2)/(3 — f) is 
integer and equal to n*. The contribution of this term 
comes from the finite s region, —Bn* ^ 
so it is not included in the scaling function f*{x). This 
term is an analytic contribution to S at f = 3— l/(n-|-l), 
n = [l..oo], including f = 5/2. 

Let us consider briefly the ordinary percolation model 
(m = 1) with an initial distribution P{s,t) decaying 
faster than s“^. Then the transition takes place at 
tc > 0, and the critical distribution P{s,tc) — 
with T = !ij2. In this situation, the Taylor expansion of 
P(s, t) around t = tc has coefficients a„ given by Eq. (1511) 
with 5/2 substituted for f. As a result, the scaling form 
of the distribution P{s,t) is 



Here we have used the fact, following from Eq. 
that the deviations from the asymptotics A„(s) — 
are of the order of si-'^+("'-i)(3-'r) fQj. 
large s. The function f*{x) is a new scaling function 
obtained from f{x) by subtracting the first n* + \ terms 
of its expansion over x’^, 

f*{x) = f{x) - ao- aix'^ - ... - a„.a;" (33) 

cr = 3 — f. In particular, if f = 5/2, then ,9 = 1, ai = 0, 
and so f*{x) = f{x) — oq. 

In Appendix iBl we find A„(s) = Bn employing the 
generating functions approach: 


P(s, t) ^ (37) 

for t approaching tc from above and below, where the 
scaling function is f{x) = oq exp(—dTruga:) in both 
phases. To obtain this result we simply replace f by 
5/2 and < by t — tc in Eqs. (1551) and (1551) . Making the 
same replacements in Eq. (1351) . we get 

^ ^ -ao E _ {^t-tc)Bx 

S 

poo 

= —|t—tc|ao J dxx~^l'^ —Ij — {t—tc)Bi 


0 if n < /3, 

= ,34) 

if n > /3, 


n -I-1 


Therefore the first sum on the right-hand side of Eq. (l32l) 
is zero, except when /3 is integer, equal to n*, and the 
sum is Bn*t^■ Thus, the singularity of the size of the 
percolation cluster is 






dxx^ ^f*{x)+Bn 


(35) 


= 47rao(|t-tc|+t-tc)- (38) 

Note that according to Eq. (IMl) . for f = 5/2, the co¬ 
efficient Bi = dxx~^/'^{e~‘^'^°‘o^ — l) = —dTrog. This 

equation describes S in both phases, giving S = 0 for 
t < tc, and S = 8ao7r(t — tc) for t > tc- 

When /3 = (r — 2)/cr is non-integer, only the scaling 
region contributes to the singularity of the percolation 
cluster size 


S = —\t — tcf / dxx^ '^/*(x). 

Jo 

In this situation, as is noted in Ref. 0, the scaling func¬ 
tion f*{x) must be different below and above tc- Namely, 
in the phase t > tc the scaling function integral 


Note that the integral J^dxx^ ^ f*{x) converges at 
the lower limit because we have subtracted from f{x) 



dxx^ '^f*{x) 


-B, 
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while in the phase t < the integral 



to comply with S' = 0. Consequently the scaling function 
f{x) = f*(x) + ao + ... + an»x^ must have a maximum 
/(a^max>0) > /(O) in the phase t < tc [1|- This asym¬ 
metry of the scaling function is observed, for example, 
for ordinary percolation at dimension d < 6 pHI ]. and 
in explosive percolation In these examples, 

0 < /3 < 1. On the other hand, when /3 is integer, the 
non-scaling region additionally contributes to the singu¬ 
larity of S. This contribution, which is important, in 
particular, for ordinary percolation above the upper crit¬ 
ical dimension 6 where /3 = 1, was not considered in 
the book [^. Thanks to the non-scaling contribution 
Bi{t — tc) to S, f(x) is the same monotonically decreas¬ 
ing function in both phases in this situation, see Eqs. dSIl) 
and (|38]) . 

B. The case of r = 3 



FIG. 2. Numerical solution of evolution equation ([3]) for 
s < 10® with the initial condition P(s, 0) = aos~^, where no = 
^(2)“^. Solid lines, curves P(s,t)s^ vs. tins for the cluster 
sizes s indicated by the numbers in the plot. Dashed line, 
scaling function from Eq. (1451) . /(tins) = ao/(l — 2aotlns). 
Inset, highlight of the small tins region, in which the curves 
here depicted approach /(tins). 


Now we show that in the case of r = 3 all derivatives 
of S{t) with respect to t are zero at t = 0 and then we 
find the respective scaling function. In this situation the 
singularity of the generating function p{z, 0) at 2 ; = I 
differs from that in Eq. 0, namely 


I - p{z, 0) = -ao(l - z) ln(I - z). 

(39) 

The inverse function is 


1 ^ 1 ^ 1 - P 

\nz -z I - 

aoln(I - p) 

(40) 

The right-hand side of this expression gives the function 
g{p) in Eq. dH]). Using Eqs. dH) and (HOI) we write the 
equation for the generating function p(z, t) near z = 1 as 
follows 

Inz = 2t(l p) + 

aoln(I - p) 

(41) 


(see Appendix [C] for the derivation). Within the radius 
of convergence of the series (HU), this expression gives the 
following scaling form of the distribution P{s,t): 

P{s,t) ^ s~^f (tins), (44) 

where 

f{x) = ao(I - 2aox)~'^. (45) 

The function f{x) plays the role of a scaling function, 
which is represented by the Taylor series (H^)) for t up to 
the radius of convergence, t <r = I/(2aolns). Surpris¬ 
ingly, in contrast to the case of f < 3, the scaling function 
for f = 3 diverges approaching Xrc = l/(2ao) from below. 
This divergence requires interpretation, since the distri¬ 
bution P{s,t) itself cannot be divergent at any cluster 
size, including s = exp[(2aot)“^]. In general, for t < r, 


At z = 1 this equation gives the relative size of the per¬ 
colation cluster for small t: 

(42) 

Therefore, if r = 3, the transition occurs at the initial 
moment, and all the derivatives of S are zero at t = 0, 
so the percolation transition is infinite-order. 

1. Scaling of P{s,t) 

In the case of r = 3, instead of the Eqs. (HHl) and m 
for f < 3, we derive the general asymptotic expression 
for the coefficients A„(s): 

(43) 


lim P{s,t)s^ ^ —>■ /(x), (46) 

s—^oo 

where x is a scaling variable such as \t—tc\s°' or \t—tc\ In s. 
If the scaling function converges everywhere in the range 
of X between 0 and infinity, then the curves P{s,t)s'^~^ 
vs. X will collapse into /(x) at sufficiently large s. In the 
case of f = 3 the curves P(s,t)s'^“^ vs. x tend to /(x) 
when s —>■ 00 in the region x < Xrc, see Fig. [5J These 
curves have a maximum near x = l/(2ao). As is shown 
in Fig. HJ the height of this maximum increases, and the 
difference between P{s,t)s‘^ and /(tins) decreases, as s 
grows. Note that the curves in Fig. [5] tend to /(tins) 
only for tins < l/(2ao), above this point they tend to 
zero. 

In the case of r = 3 we cannot formulate a relation 
similar to Eq. dSlD between the scaling function /(x) and 


A„(s) ^ 2"(n+l)a”+is-2 (Ins)" 










the singularity of S{t) at 0 . All the coefficients Bn of 
the series S{t) = — X]n>i 0, which indicates an 

infinitely smooth singularity of S{t) at 0. 


C. The singularity of susceptibility 

Let us find the singularity of the susceptibility and its 
relation with the percolation cluster size S. In mean- 
field models, near the critical point, the susceptibility 
typically follows the Curie-Weiss law y ^ |t — It 

is easy to see that this is the case for the ordinary per¬ 
colation model if the exponent of the initial cluster size 
distribution f 7^ 3. In ordinary percolation the suscep¬ 
tibility is the average cluster size to which a uniformly 
randomly chosen node belongs 0, namely 

X = (s)p = (47) 

S 

For m = 1 the singularities of y and S can be related in 
the following way. Summing over s both sides of Eq. m 
gives 


'm 


2S{s)p = 2Sx- 


(48) 


In the case of r < 3 the percolation cluster size S has a 
power-law singularity with exponent (3 = (f — 2)/(3 — r), 
Eq. (fTTI) . then 


\d\nS ^ t-2 

2 dt ~ 2(3-f) 


(49) 


for small t. 

At f = 3, we obtain an anomalous singularity of the 
susceptibility which diverges at t = 0, y(t=0) ~ 
Substituting S ^ g-i/(2aot) into Eq. (H5)) we 

find 


l oins' _ J_ 
2 dt 4ao 


(50) 


differing from the Curie-Weiss law. Note that when 
7 > 3, the initial susceptibility does not diverge, and 
the transition occurs at tc > 0 with standard exponents. 
When f < 3, the initial susceptibility diverges, the transi¬ 
tion occurs at tc = 0, and the susceptibility demonstrates 
the standard behavior, y ^ 


IV. EFFECT OF INITIAL CONDITIONS ON 
EXPLOSIVE PERCOLATION (m > 1) 

In this section we extend the analysis made in the pre¬ 
vious section for m = 1 to to > 1. We consider the 
following explosive percolation model BE mi- The 
evolution starts from a given distribution of clusters. At 
each step, we choose at random two set of to nodes. Then 
the node in the smallest cluster of each set is selected. 


and these two nodes are interconnected. The evolution 
equation for an arbitrary to takes the form: 


dP{s, t) 


s-l 

= s ^ Q{u, t)Q{s-u, t) - 2sQ{s, t), 

U — 1 


(51) 


where Q{s,t) is the probability that a chosen node seats 
in a cluster of size s. The relation between Q{s) and P{s) 
is 


Q{s) = 


s-l 


i-Y.p{u) 


u—1 


l-^P(n) 


= mP{s) 


U—1 
-\ m—1 


l-^P(n) 


(52) 


where the last approximate equality takes place at large 
s B]. In the phase with the percolation cluster, this 
relation between distributions P(s,t) and Q{s,t) is sim¬ 
plified: 


Q{s) ^ TOS'™-ip(s). (53) 


Then we can write the partial differential equation for 
the generating function p{z, t) defined in Eq. (|3]) for any 
to: 


= -2[5(t)]2(—1) K(l-p) - m{m-l)S{t)\ 

"(54) 

This partial differential equation can be solved by ap¬ 
plying the hodograph transformation 2l| ) which leads to 
the following equation 

/■* , 91nz 

1= 


= / dt'2[S(t')Y^'^ [to^(1 —p) — to(to— l)5'(t')] 


+5(p), 


(55) 


where the function g{p) is determined by initial condi¬ 
tions. Similarly to Eq. ®, we have 


9{p) = - - 


1 -p 


l/(f-2) 


(56) 


aor(2 - f) 

At 0 = 1, the generating function p(l,t) = 1 — S[t), so 


0 = 2m'^S{t) [ 

Jo 


- 2 ( to "- to ) f dt'[S{t')]'^"'-^ 
Jo 


S{t) 


1 i/f-f-a) 


aor(2-f) 


One can see that 


S = Btf^ 


(57) 

(58) 



















is a non-trivial solution of Eq. ((Ffl) with 
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and 


/3 


t-2 

1 — (2m — l)(f — 2) ’ 


(59) 


way to the case of m = 1. Using Eq. (EU at t = 0 we 
write 


s-l 

Ai(s) =s ^ Q{u, Q)Q{s-u, 0) - 2sQ{s, 0) 

U—1 


B = [-aor(2-f)] 


/3/U-2) 


2m 


(f-2)[l+(m-l)(f-2)]' 


(3 - f)/3 


(60) 

In Ref. we introduced the generalized susceptibility 
for these explosive percolation models. The susceptibil¬ 
ity X is introduced in terms of the probability C2 that 
two nodes selected by our algorithm belong to the same 
cluster 


C2 = 


^2m 


N 


J_ sQ'^{s) 2m 
P{s) ^ 


(61) 


s-l 

=S “ '3(s, 0)] [Q(s-u, 0) - Q{s, 0)] 

U—1 


oo 

-s{s - l)Q{s, 0)2 - 2sQ(s, 0) Q{u, 0), (62) 

U—S 


where we have used the normalization condition 
Q(s,0) = 1. The asymptotics of Q(s,0) are obtained 
substituting the power-law initial distribution P(s, 0) = 
Ao(s) = aoS^~^ into Eq. (15^ : 


In particular, when m = 1, this susceptibility is reduced 
to the standard one for ordinary percolation y = (s) p [^ . 
In [HI we showed that the susceptibility is divergent at 
t = 0 if the exponent of the initial cluster size distribution 
T < 2-1- l/(2m— 1). For f > 2-|- 1/(2 to — 1) the transition 
occurs at tc > 0, with the critical exponents and scaling 
functions calculated in P,[2l|. If b < 2 + l/{2m- 1) 
the size of the percolation cluster follows the power-law 
S = with /3 and B given by Eqs. (IMl) and (1501) 
respectively. Here the case of f = 2 -|- l/(2m — 1) for 
m > 1 distinguishes itself significantly from that for m = 
1. In the next subsection we show that an infinite-order 
percolation transition takes place at tc exceeding zero 
when f = 2 -I- l/(2m — 1) and m > 1. 

The remainder of this section is organized in the fol¬ 
lowing way. In Sec. IIV Al we consider the region t < 
2 -I- l/(2m — 1). For this range of r we derive the scal¬ 
ing form of the distribution P(s, t) and relate the scaling 
function and the singularity of 5. In Sec. IIVBI we con¬ 
sider the case of f = 2 -I- l/(2m — I) and show that in 
this situation we have an infinite-order phase transition. 
Next we derive scaling functions for this problem and the 
singularity of the relative size of the percolation cluster 
S. Finally in Sec. IIV O we find the critical singularity of 
susceptibility. 

A. The case of r < 2 -|- ll{2ni — 1) 

Let us first consider the region f < 2 -f l/(2m — 1), 
derive the scaling form of the distribution P{s,t), and 
obtain the critical singularity of the percolation cluster 
size. 


Q(s, 0) = mAo{s) 


Y^o{u) 


'l^S 


m— 1 


Then the asymptotics of Ai is 

2m(2-f) 


m(2 —f) —1 
mf_ 

° (f-2)™-!’ 

(63) 


2^2m 


Ai(s) =m^aQ 


-l-2y 


(^_ 2 ) 2 m -2 

+ [ dx[x™(2-U-i_i][(i_a.)-(2-U-i_i] 

^0 


=ais2™(2-U 


where 


Oi = 


ao™m^r[—m(T — 2)]^ 
(f_2)2m-2p[_2^(f _ 2)]' 


(64) 


(65) 


For f = 2 -I- l/(2m) the coefficient oi becomes zero, and 
we need to find the next term. This term can be obtained 
by taking into account the next-to-leading order term in 
the expansion of the right-hand side of Eq. (EH) in powers 
of s. For this sake we use Eq. (fTBl) . In this way, when 
Oi = 0, we get the following asymptotics of Ai(s) oc 

gm(2-f)-l ^ g-3/2 

Repeating the procedure for the next coefficients 2l„(s) 
we find that the n-th coefhcient has the power-law 
asymptotics: 


A„(s) ^ (66) 

where the exponent 


1. Scaling of P{s,t) 


(7 = 1 — (2m — 1 )(t — 2), (67) 


For m > 1, the coefficients A„(s) of the Taylor expan¬ 
sion of P{s,t), Eq. (IT^ . can be calculated in a similar 


for any m. We express the prefactors a„ in terms of m, 
r and oq, similar to ai. For instance, for 02 and 03 we 
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find: 

2a^"^~^m^r[l-m(f-2)]3r[2-(3m-l)(f-2)] 

“2 (f _2)4m-2r[2-2m(f-2)]r[l-(4m-l)(f-2)]’ ^ ’ 


4ag™“Wr[l-m(f-2)]4r[3-(5m-2)(f-2)] 

3(^_2)6m-4r[2_(6m-2)(f-2)]r[2-2TO(f-2)]2 


X 


r[2-(3m- l)(f-2)] 


r[2-(3m- l)(f-2)] 
r[3-(5m- 2)(f-2)] 


[l-(4TO-l)(f-2)]r[2-2TO(f-2)] \ 
TO(f — 2)r[3— (4m— l)(f—2)] j 


— (m—l)(f— 2)r[— m(f— 2)] 


(69) 


When m > 1, we have to derive expressions for the co¬ 
efficients a„ individually, unlike the general expression 
(EU for all a„ in the case of m = 1. Each coefficient 
an>o becomes zero when f = 2 + n/[n(2m—1)-|-1]. In 
this case, the resulting asymptotics of ^„(s) decay as 

g-l-[n(m-l) + l]/[n(2m-l) + l] 

The sums the equations for 

^ri>2(s), are obtained by the generating function ap¬ 
proach explained in Appendix |B] in detail for m = 1. 
We find the singularity 

Pn{z) = a„r[2-f-|-ntT](l-2:)'^“^“”‘^, (70) 


by differentiating the evolution equation m and rela¬ 
tion (l5^ . The value of = /9„(1) diverges or converges 
depending on f, similarly to the case of m = 1: 


B 


n 


0 if n < /3 

a„r[2— f-|-n{l—(2m—l)(f— 2)}] if n = /3 


oo 


if n > /3 


(71) 

where /3 = {f — 2)/[l — (2m — l)(f — 2)]. 

The function P(s, t) can be written in the scaling form: 


Pis, S) ^ s^-^ Y. = s'-V(s<5'/"). (72) 

n=0 


Here we used the asymptotic behavior of the coefficients 
A„(s). In the vicinity of x = 0, the expansion of the 
function /(x) is 

/(x) = Oo -|- oix”^ -I- a2X^‘^ -I- aax^'^..., (73) 


where a = 1—(2m—l)(f—2) and the coefficients a„ are 
given by Eqs. (1551) . (15^ . (15^ . etc. The scaling function 
and S are interrelated in the same way as in Eq. (IS3 
from Sec. inEi 


dxx^ '^f*ix)+B„ 


(74) 


where cr = 1 — (2m — l)(f — 2), /3 = (f — 2)/cr, and 
n* = . The function /*(x) = /(x) - Onx""^ as 

in Eq. (l33l) . The constant 

Bn- =^A„.(s) =p„.(l) = a„.r (2-T 4-71*0-) 

S 

in nonzero only when /3 is integer, /3 = [/3J = n*, see 
Eq. (1711) . 


B. The case of r = 2 4-l/(2m — 1) 

Let us consider the case of f = 2-|-l/(2m—1) and show 
that in this situation we have an infinite-order phase tran¬ 
sition. When m > 1, the coefficients A„ in the case of 
f = 24-1/ (2m — 1) have exactly the same form (l66l) as for 
f < 2-|-l/(2m—1). This coincidence is due to the conver¬ 
gent convolution integral in Eq. (IM)) for Ai and similar 
integrals for A„>i. This is in contrast to m = 1, f = 3, 
where we have the coefficients An ^ s~^ (Ins)" with log¬ 
arithmic factors unlike An Jqj. rn = 1, f < 3, 

see Sec. IIIIBI Recall that these logarithms emerged due 
to the divergent convolution integral in Eq. (1151) for Ai 
at T = 3 and similar integrals for A„>i. In this respect, 
there is a principle difference between m = 1 and m > 1. 
Below we show that for m > 1 and f = 2 4- l/(2m — 1) 
an infinite-order percolation transition takes place at a 
tc > 0 . 


1. Expansion of P{s,t) att = 0 

When -7 = 24-1/(2m — 1), we have the exponent cr = 0, 
so the coefficients A„(s) oc s^”"^ for all n, see Eq. (1661) . 
Consequently the asymptotics of the cluster size distri¬ 
bution is 

OO 

Pis, t) ^ s-2W(2m-l) ^ ^-2m/{2m-l) j 

n—0 

(75) 

Substituting T —2 = l/(2m—1) into Eqs. (1551) . (|55)) . (1551) . 
etc., we find the general expression for the coefficients a„ 

aor[n-|-l/(2m—1)] 
n!r[l/(2m-l)] 

2m[ao(2m—l)]^™“^r[(m—l)/(2m—1)]^ 
r[2(m-l)/(2m-l)] 

(76) 

and so 

/(t) = ^a„t"=ao(l--) , (77) 

where 

(2m-l)i-2™r[2(m-l)/(2m-l)] 
2a^'"-imr[(m-l)/(2m-l)]2 






(78) 
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FIG. 3. Numerical solution of the evolution equation dsa, 
for m — 2 and s < 10® with the initial condition P{s, 0) = 
aos~^^^, where uq = ^(4/3)“^. Solid lines, curves P{s,t)s‘^^^ 
vs. t for the cluster sizes s indicated by the numbers in the 
plot. Dashed line, function f{t) in Eq. (1771) . 


The function f{t) diverges at t = r. Fig. [3] demonstrates 
how the curves P(s, vs. t approach f{t) in 

the region t < r as s approaches infinity. 

Substituting Eq. (ESI) into Eq. (ESI), we obtain the 
asymptotics of the distribution Q(s,t)\ 

Q{s,t) ^ TO( 2 m-l)™-V(f)'"s-i-'"/( 2 m-i)^ 

The first moment of this distribution 

(s)q = ^ sQ{s, t) ^ f{tr E 

S S 

is divergent in the interval 0 < t < r. Note that also 
(s)p = f{t) Es ^ ^-1 gj.g 

divergent in this interval. Summing over s both sides of 
Eq. (|5T1 i we get 


^ = 25-(s)q. (79) 

Due to the divergence of {s)q, the only solution of this 
equation in the interval t < r is S' = 0. Thus we have a 
phase without percolation and with divergent susceptibil¬ 
ity, which enables us to call this phase “critical”. Below 
we show that the percolation threshold of this transi¬ 
tion is exactly r. For t > r, we show in Appendix iFl that 
the moments {s)q and (s)p are finite. 


2. Scaling of P{s,t) at tc 

For t <tc the function f{t) is convergent, and for large 
s the curves P(s, collapse into the func¬ 

tion f{t), see Fig. |3 The value of P(s, 
approaches infinity as s —>■ oo. Our numerical results. 
Fig. m indicate that P(s, grows linearly 



FIG. 4. Logarithmic contribution to the asymptotics 
of P{s, tc) for (a) m = 2 and (b) m = 3. Solid lines, 
P(s, for s from 1 to 10® found numerically from 

the evolution equations with the initial condition P(s, 0) = 
where ao = ^[2m/(2m — 1)]”^. Dashed lines, 
straight lines presented for reference. 


with Ins for sufficiently large s. In the range s S [1..10®], 
however, it may be difficult to distinguish different slowly 
varying functions, such as powers of logarithm. There¬ 
fore, we cannot exclude the possibility that the solid line 
asymptotics follows a law (Ins)^ with exponent A close 
to but different from 1. 

In Appendix [D] we obtain the following scaling form of 
P(s, t) near 1^. 


P{s,t) ^ C(lns)^s-2'"/(2'"-i)/ [C'(lns)^(2m-i)(^_^j 


(80) 

where C is a constant, A > l/(2m—1), and f{x) is the 
function defined in Eq. EH)- Then, the critical distribu¬ 
tion P{s,t = tc) behaves as 


P(s, tc) = aoC'(ln s)As-2m/(2m-l) 
for large s. Recall that P(s,0) = 
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3. Critical behavior of S 


In Appendix |E] we derive the critical singularity of the 
percolation cluster size S: 


S ^ exp(—dd 
where d = t — tc > 0, 

1 


d = 


A(2m — 1) — 1 ’ 


(82) 


(83) 


and d is a constant 



/m(3m—2)[A—l/(2m—1)]^ 
y (m—1){—ooCT[—1/(2 to—1)]}^“^™ 


1/[1-A(2m-1)] 


(84) 


Note that A cannot be smaller than l/(2m — 1). 


C. The singularity of susceptibility 

Let us we find the critical singularity of susceptibility. 
The critical behavior of y is determined by large s. In 
the region in which S' > 0, i.e. t > tc, we can substitute 
(5(s) = mS"^~^P{s), Eq. (IMl) . into Eq. (IbTI) . which gives 

y ^ sP{s, t) = m^S^^-\s)p 

S 

close to tc- Summing Eq. (EU over s gives 
oc 

_ = 25”^ y sQ{s, t) ^ 2mS‘^^-^{s)p. 
at 

s 

Combining the last two equations we finally get 

^ m dS md\nS 
^ ^ ^ T dt ■ 

This remarkably general formula relates the critical sin¬ 
gularities of the susceptibility and the percolation cluster 
size in all situations shown in TablelH These situations in¬ 
clude the finite- and infinite-order continuous phase tran¬ 
sitions. 

For f < 2 + l/(2m—1), when S = Bt^, this equation 
ensures that the susceptibility has the following singular¬ 
ity: 


(85) 


( 86 ) 


(87) 



( 88 ) 


For m > 1 and f = 2 -|- l/(2m—1) the susceptibility 
is divergent below while above tc it has a power-law 
singularity 


X = 




(89) 


which we obtained by substituting Eq. (15^ into Eq. (EZl). 



FIG. 5. Critical behavior of the susceptibility y for (a) m = 
1, (b) m = 2, and (c) m = 3. Solid lines, numerical solution 
of 10® evolution equations with initial condition P{s, 0) = 
(^[2m/(2m—Dashed lines, (a) power-law 
in Eq. (1501) : (b) and (c) power-law in Eq. (18911 assuming p = 
l/(2m — 2), and C = 0.024 and 0.0088 for m = 2 and 3, 
respectively. 


Figure [5] presents the susceptibility of the infinite-order 
percolation transitions for f = 2 -|- l/(2m—1). We ob¬ 
tain the solid curves in this figure inserting the numeri¬ 
cal solution of 10® evolution equations into the expression 
y = sQ(s, tY/P{s, t). The dashed lines are the power 
laws in Eq. o and Eq. (IMl) for m = 1 and m > 1, 
respectively. For m = 2 and 3 we used the value of 
/r = l/(2m — 2) corresponding to the asymptotic distri- 
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bution P{s,tc) s 2m/(2m 1) which is observed in 

Fig. m 


V. CONCLUSIONS 

In this article we have explored the impact of the initial 
cluster size distribution in a set of models generalizing or¬ 
dinary percolation. Specifically, we focused on explosive 
percolation models, but our approach could be applied 
to a much wider range of generalized percolation mod¬ 
els that can be reduced to various aggregation processes. 
In particular, we considered initial cluster size distribu¬ 
tions, for which the percolation phase transition turned 
out to be remarkably different from that for the evolution 
started from isolated nodes [l3, [H, [2l| . Our results are 
summarized in Table ID We have found the special values 
of the exponent f of the initial cluster size distribution 
s~'^, namely f = 2 + I/(2m — 1), for which the perco¬ 
lation cluster emerges continuously, with all derivatives 
zero. For ordinary percolation, m = I, this singularity is 
at t = 0, and the susceptibility diverges as in contrast 
to the Curie-Weiss law typical for mean-field theories in¬ 
cluding various percolation problems, in which the evo¬ 
lution starts from isolated nodes. We have found that for 
explosive percolation, m > 1, the situation is even more 
interesting. When f = 2 + 1/{2m — 1), (i) the phase tran¬ 
sition occurs at tc > 0; (ii) the transition is continuous, 
of infinite order; (hi) the phase 0 < t < tc is critical in 
the sense that the generalized susceptibility diverges at 
any t < and the size distribution of clusters has the 
same asymptotics s~'^ in the entire critical phase; (iv) the 
susceptibility diverges above the transition with a criti¬ 
cal exponent different from 1; (v) the size distribution of 
clusters at tc is the same power law but with additional 
logarithmic factor, see TableHl Finally, in the special case 
of r = 2 -I- l/(2m — 1), we have obtained unusual scaling 
both for ordinary and explosive percolation. Note that, 
counterintuitively, in the case of explosive percolation, 
for this special value of t the power-law critical singular¬ 
ity of the generalized susceptibility is accompanied by a 
strong divergence of the moments (s)p and (s)q at the 
critical point, (s)p ^ (s)q ^ exp[const(t — tc)~^], see 
Appendix [F] We studied the range f > 2, in which the 
initial distribution P(s, 0) ~ 3^“"^ is normalizable and 
the average cluster size is finite. In the case of f < 2 and 
TO = 1 it was found that S{t > 0) = I [^ . 

The infinite order singularity for the per colation clus¬ 
ter at t = 0 was found in Ref. 1^. 1^ i n aggrega¬ 
tion processes with power-law kernels [2 instead of 
our power-law initial distribution of clusters for ordinary 
percolation. This singularity at zero was also observed 
in epidemic models and percolation problems on equi¬ 
librium scale-free networks with the degree distribution 

P((7)^g-3 iii,[2|. 

The critical phase and the inhnite-order phase transi¬ 
tion resemble the Berezinskii-Kosterlitz-Thouless transi¬ 
tion [sO, [Ml , which was observed in numerous systems 


at a lower critical dimension. In addition, these singu¬ 
larities were observed in heterogeneous one-dimensional 
systems with long-range interactions [s^, , in various 

growing networks [3J-l39|, and in percolation on hierar¬ 
chical and nonamenable graphs Interestingly, in 

these growing networks, the critical distributions of finite 
cluster sizes also had factors with powers of logarithms 


In conclusion, we have found that initial conditions 
can have dramatic effect on the scenario of percolation 
transition and its various generalizations. For particu¬ 
lar initial distributions of clusters in the case of explo¬ 
sive percolation, we revealed a continuous phase tran¬ 
sition of infinite-order singularity, which resembles the 
Berezinskii-Kosterlitz-Thouless transition. We suggest 
that our findings are valid for a wide range of gener¬ 
alizations of percolation, in particular, for explosive per¬ 
colation models with power-law kernels. 
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Appendix A: Calculation of the coefficient A 2 (s) for 
m = 1 and r < 3 

Let us derive expression mi) for the coefficient ^2(5) 
of the expansion m of P{s,t) in the case of ordinary 
percolation. We differentiate both sides of Eq. m with 
respect to t and replace P{s,t)\t=o with Ai(s)i!. In 
this way we get 


s-l 

^2(5) = s ^ Ai(m)Ao(s— u) — sAi(s). (Al) 

U—1 

Rearranging this equation in order to cancel divergencies, 
we write 

s-l 

^2(5) = s ^ - Ai(s)] [Ao(s-u) - Ao(s)] 

U—1 

CO 

- s(s - l)Ao(s)Ai(s) - sAi(s) ^ Ao(u) 

U—S 

S-l 

+ 'sAo(s) Ai{u). (-A.2) 

Li —1 

Similarly to Ai, let us replace the sums over u in Eq. (IA2I) 
by convergent integrals. The first and second sums on the 
right-hand site of this relation can be directly replaced 
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by the respective integrals. So we have to analyze only 
the third sum. For this sum, there are three possibilities: 
(i) if f < 5/2, then the asymptotics of Ai{u) have the 
exponent 4 — 2-? > —1, (ii) if f > 5/2, then the expo¬ 
nent 4 — 2f < —1, and (hi) if f = 5/2, then the coeffi¬ 
cient oi from Eq. m is zero, and Ai(s) decays as s 
see Eq. (Ha- In the the first case the sum 
can be directly replaced by an integral convergent at the 
lower limit. In the second and third cases this sum must 
be first replaced with ~ Su>s^i(^)' The 

sum is a finite constant and J2u>s ^i('a) can 

be replaced by a convergent integral. Substituting the 
asymptotics of Aq and Ai , we obtain 


^2(5) = aoais^"^" [ dx - l] [(l-a;)^“" - l] 

^0 


ooUis^ — uoais"^ 


dx c 


+ 0(5/2—f)aoais^ / dxx'^ 


+ 9(f—5/2)aos‘^ 


—ois^ dxx'^ 
.u=l dl 


(A3) 


where the step function 9(x) is defined here as 0(x<O)=O 
and 9(x>0)—l. When f>5/2, the fourth term on the 
right-hand side of this equation is zero, and in the last 
term, the sum '^^Ai{u) = 0, see Appendix [BJ When 
f=5/2, the coefficient ai=0, and so the asymptotics of 
^2(5) is given '^c show in Ap¬ 

pendix [B] in this special case the sum J2u ^i('a) is finite, 
see Eq. (lB7l) . As a result, for any f < 3, we have 


tion (m by and summing over s, 

's-l 


Piiz)=2Y^ 

S 


z^uAq(u)z'^ “Ao(s—m) 

_U—1 


0'*sAo(s) 


= 2 0'’sAo(s) 

S 


^z“Ao(u)-l 


= dim{po(z) - 1)^. 


(B2) 


Here we used Xls fiu)f(s-u) = J2s /(«) E« f{u). 

The function po(z) is the generating function of the 
initial distribution P(s,0) = Ao(s). For a power-law 
P(s, 0) = the singular behavior of Pq(z) at z = 1 

is given by Eq. which we reproduce here for the sake 
of clarity. 


l-po(z) ^ -aor(2-h)(-lnz)"-2. (B3) 

Inserting this result into Eq. (IB2I) we find the singularity 
of pi(z) at z = 1: 

pi(z) ^ -2alT(2 - ff{f - 2)(1 - zf^-\ (B4) 

Differentiating both sides of the evolution equation with 
respect to t, and combining with Eq. (IB2I) and (IB3I) . we 
express the function p 2 (z) as 

P2(z) =Y^ z‘'Ao(s) y] z“mAi(u) 

S U 

+Y.^-ms)Y. z^uAoiu) - y]z"sAi(s) 

s u s 

= ^d\r,z [d\nz(po(z) - 1)^] 

^2a^r(2 - Tf{f - 2)(3f - 7)(1 - z)^^-®. (B5) 


A2(s) 


3 2F(2-f)3 
3F(6-3f) 


s 


7-3f 


(A4) 


Appendix B: Generating functions approach for 
ordinary percolation (m = 1) 


Using generating functions, we will obtain 
for different exponents f and the general expression for 
the coefficients a„. Let us consider the generating func¬ 
tions of the coefficients A„(s) in the series (HU, that is 


Pn(z) =^z'‘An(s). (BI) 

s 


We obtain the generating function pi(z) in terms 
of po(z) multiplying both sides of the evolution equa- 


In a similar way, for a general n, we obtain 

_2"a”+ir(2 - f)'‘+ir[(n -f I)(3 - h) - I] 
(n+I)!r[(n + I)(2-f)] 

x(I-z)(”+i)(^-®)+^ (B6) 

The sums ^n(®) equal to the value of the 

generating function p„(z) at z = I. Then according to 
the last equation 


B 


n 


0 


2"[-aoF(2-f)]’"+i 
n -I- I 
oo 


if n < /3, 


if n = /3, 
if n > /3, 


(B7) 


where /3 = (f — 2)/(3 — t) is defined by the condition 
1 — T -I- /3(3 — t) = —I. 

Equation (IB6I) can be used to find the general form 
of the prefactors a„ in the asymptotics of A„(s). We 
obtain the singular contribution to Pn(z) near z = I by 
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inserting the respective power-law An(s) = r+na Appendix D: Expansion of P{s,t) at tc > 0 for m > 1 

into Eg. (IBII) : and f = 2-I-l/(2m-l) 


Pn{z) ^ a„r[(n+ 1)(3-f) - 1](1 - z)(-+i)(--3)+i. 

Comparing Eqs. (jB8l) and (IB6I1 we readily get: 

[2aor(2 - f)]"+' 
2r[(n+l)(2-f)](n + l)!- 


(B8) 


(B9) 


Appendix C: Calculation of coefficients A„(s) for 
m = 1 and f = 3 


For the particular case of m = 1 and f = 3 we will cal¬ 
culate the asymptotics of the coefficients A„(s), which 
differ from those in Sec. IIII Al for f < 3. The difference 
follows from the divergence of the convolution integrals 
of Eqs. (ITSl) . (IA3I) . etc, for the coefficients Ai(s), ^2(5), 
etc., when f = 3. The arguments of these integrals di¬ 
verge at the upper and lower limits as (1—a;)“^ and 
respectively. To remove these divergences we must sub¬ 
tract one extra term to each factor of the convolution 
argument. For Ai, instead of Eq. (fHl) and (ITSl) . we write 

s-l 

24i(s) = s [Ao(it) — aos~^ — 2aos~^(s — u)] 

U—1 

X [Ao(s — u) — aoS~^ — 2aoS“^u] 

s-l 

— dogS”® ^ u{s — u) — ao(s — l)s“^ 

U—1 

S-l s-l 

— dOgS”^ ^ U + dags“^ ^ uAq{u) 

u—1 U^l 


S-l 


-I- 2agS ^ ^ Ag(u) — 2sAg(s) 


u,= l 


= ttnS 


' f dx{ [x ^ — 1 — 2(1 — x)] 
Jo 


X [(1 — x) ^ — 1 — 2x] — dx(l — x) — dx} — OgS ^ 

/•OO 1 

dxx“^ -I- dogs”^ ^ (Cl) 


— 2anS ^ 


Since X]«<s “ ^ ^ term decays as s ^Ins, 

while the others decay as s“^. Then, for large enough s, 
we obtain 


Ai(s) = dcgS ^Ins. (C2) 

We obtain the asymptotics of the other coefficients A„(s) 
in a similar way. The An found by this procedure have 
the general expression 

An{s) ^ 2"(n+l)a((+is-2 (In sf . (C3) 


Let us obtain the scaling form of P{s,t) near the 
critical point tc > 0 in the case of m > 1 and f = 

2 -I- l/(2m—1). Our numerical solution of exact evolu¬ 
tion equations shows that that asymptotics of the distri¬ 
bution P{s, tc) differs from a pure power-law by a factor 
In^ s with A close to 1, see Fig.|3]and Appendix lEl Then, 
at tc we have 

P{s,tc) = agC'(lns)^s-2'"/(2'"-i)^ (pi) 

where C is a constant. 

By the same method as in the Secs. IIII A[ IIIIBIIIV Al 
and II V B1 we obtain the expansion of P{s,t) around tc 
starting from the distribution P{s,tc) in Eq. (IDll) . The 
resulting asymptotics of the coefficient A„ is 

A„(s) ^ a„[C(lns)^]i+”(2™-i)s-2W(2m-i) ^ 2 ) 

with prefactors a„ given by Eq. Gil). The series with 
these coefficients gives 

P{s,t) = C(lns)^s- 2 '"/( 2 '«-i)j [(i _ t^){C(lns)^}2'"-i] , 

(D3) 

where /(x) is the function from Eq. (177)) . This scaling 
form is similar to Eq. (jddll for m = 1, f = 3, that is, 
a logarithmic factor of s appears in the argument of the 
scaling function. The difference, however, is that for m = 

1 we have tc = 0, while for m > 1 we have tc > 0. 


Appendix E: Percolation cluster size singularity at 

tc > 0 for m > 1 and f = 2 -|- l/(2m—1) 


In this appendix, we show for m > 1 and r = 
2-|-l/(2m—I) that an infinite-order percolation transition 
takes place at the critical point tc = r given by Eq. dZH]). 
When t < tc, the percolation cluster is absent. S' = 0. 
Let us obtain the critical singularity of S(t > tc) and 
relate it with the asymptotics of P{s,tc). 

Figure m shows that at tc the asymptotics of P{s,tc) 
differs from the initial power-law one by a factor 


H{s) 


P{s,tc) 

aos-2”i/(2m-i)' 


(El) 


For the sake of brevity we do not indicate that the 
function H{s) depends on m. Substituting P[s,tc) = 
iL(s)ags“^™/(^™“^) into Eq. dT]) (this equation is valid 
for any m) we find the asymptotic function near z = 1 
at the critical point, i.e., p{z,tc) = P{s,tc)z‘’ 

1 - p{z, 4) ^ ag ^(1 - Z*)iL(s)s- 2 W( 2 m-l) 


^ ag(-lnz)l/(2™-l)jr 

^ -aor[-l/(2m - 1)]H (- Inz)^/^"'"-^) . (E2) 
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Here we assumed that H(s) grows with s 
slower than any power law, which means that 
H{—x/ \nz)/H{—l/\nz) —>■ 1 when z —?■ 1. 

The inverse function of the singularity of p{z,tc) is 
the the initial condition for solving the differential equa¬ 
tion (l54|) in the phase with S' > 0. Inverting Eq. (|E2ll we 
get 


( _ ^ _ 

Vaor[-l/(2m - 1)]H [(1 - 


2m —1 


(E3) 


at t = tc- With this initial condition, the solution of 
Eq. (IMl) at z = 1 is 


_ zM _ 

aor[-l/(2m-l)]iS [S(t)i-2'"] 


2m—1 


2m‘^S{t)f (E4) 

Jtc Jtc 

for t approaching tc from above. 

Let us introduce the function 


for — InS = H. With the initial condition S{tc) = 0 the 
solution of this equation is 


S ~ exp(—(ii5 


where 


1 

A(2m- 1) - 1’ 


and d is a constant 


(ElO) 


(Ell) 


, 1 f w(3m-2)(m-l) ^ \ ^ 

2m-l\p{-aoCr[-l/{2m-l)]y-^^ J ^ 


The values of p from 0 to oo corresponds to exponent 
A, respectively, from oo to l/(2m—1). The latter is the 
lower bound for A. 


Appendix F: First moments of distributions P{s, t) 
and Q{s, t) 


D{t) = -\nS{t). (E5) 


For small t — tc the integrals in Eq. (IE4I) can be calculated 
by changing the integration variable to nD(t') = x, where 
n = 2m—2 and 2 to—1 for the first and the second integral 
on the right-hand side of Eq. (IE4|) , respectively. 



dt' exp[—nD(t')] 


D(t) 


dx 



-1 

exp(—x) 




nD{i) 


dx exp(—x) 




(E6) 


Notice that the lower limit of the integral in this equation 
nD{t) —^ oo when t approaches tc- Substituting Eq. (IE6I) 
into Eq. (lE4l) we obtain 


H ^ _ (2m-l)(m-l) dD{t) 

{—floT)—l/(2m—l)]}i“2'" m{3m — 2) dt 

,(E7) 

Our numerical solution of exact evolution equations, 
Fig. a shows that at the critical point, P(s,tc) 
s^”'^ln'^s, where the exponent A was found to be close 
to 1. This asymptotics implies that 


H{x)^C{\Yixf, (E8) 


for large x, where C is a constant. Inserting this expres¬ 
sion into Eq. (lE7l) we obtain the differential equation 

dP m(3m-2){-aoCr[-l/(2m-l)]}^—1 
dt (to— l)(2m—’ 

(E9) 


In this section we find the singularities of the first mo¬ 
ments of the distributions P(s,t) and Q{s^t), {s)p and 
(s)q, respectively, using the results of the previous sec¬ 
tions. For TO = 1 the distribution Q{s,t) = P{s,t), and 
so (s)p = {s)q = X- In the region f < 3, according to 
Eq. (H^. we have 

2(3-fPt 

(Fl) 

while for f = 3 Eq. (1501) gives 



(F2) 

For TO > 1 we find the singularities of (s)p and {s)q 
inserting the previously obtained singularity of the rel¬ 
ative size of the percolation cluster S into Eqs. (l86l) 
and (1701). respectively. For f < 2-|-1/(2 to— 1) the power 
law S = Bt^, Eq. (1551). results in the first moments 

d2 (1—m) 

{s)p ~ ^ i2(l-m)/3-l 

2m 

(F3) 

and 


d1 —m 

(5)q = 

(F4) 


where (3 and B are given by Eqs. dsni) and O, respec¬ 
tively. For TO > 1 and f = 2-|-1/(2 to— 1), the infinite- 
order transition occurs at tc > 0 with the singularity 
S ^ exp{—dS~^). This critical behavior leads to the fob 


lowing singularities of the first moments: 



(s)p - 

j exp[2d(TO — 1)<5 

(F5) 

and 





{s)q ' 

^ exp[d{m — 

(F6) 
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for t > tc, where the exponent /i and the constant d are 
given by Eqs. (1551) and (1551) . respectively. In this situa¬ 
tion the first moments (s)p and (s)q exhibit exponential 


divergences approaching tc from above. In the “critical 
phase”, t < tc, the moments (s)p and (s)q and the sus¬ 
ceptibility X are divergent, see Sec. lIVBl 
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